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ABSTRACT 

We present the continuation of our long-term spectroscopic monitoring of the gravitationally lensed quasar QSO 2237-F0305. We 
investigate the chromatic variations observed in the UV/optical continuum of both quasar images A and B, and compare them with 
numerical simulations to infer the energy profile of the quasar accretion disk. Our procedure combines the microlensing ray-shooting 
technique with Bayesian analysis, and derives probability distributions for the source sizes as a function of wavelength. We find that 
the effective caustic crossing timescale is 4.0 +1.0 months. Using a robust prior on the effective transverse velocity, we find that the 
source responsible for the UV/optical continuum has an energy profile well reproduced by a power-law R oc with f = 1.2 ± 0.3, 
where R is the source size responsible for the emission at wavelength A. This is the first accurate, model-independent determination 
of the energy profile of a quasar accretion disk on such small scales. 

Key words. Gravitational lensing: quasar, microlensing — Quasars: general. Quasars: individual QSO 2237+0305, Einstein Cross 



1. Introduction 

Soon after quasars were discovered, it was suggested that they 
are powered by the accretion of gas on to supermassive black 
holes at the centers of galaxies. Since the infalling matter has 
non-zero angular momentum, it forms a disk spinning around the 
central black hole (Lvnden-Bell 1969). Viscosity within the disk 
is thought to result in an outward transfer of angular momentum, 
thereby allowing the material to spiral gradually inward, heating 
the disk and causing it to emit intensely over a wide range of 
wavelengths (IShieldslll9'78l) . 

Despite the numerous studies addressing the subject, there 
are still very few direct observational constraints on the spa- 
tial structure of quasar accretion disks. Quasars are located at 
cosmological distances, and it is particularly difficult to observe 
the inner regions of these objects. Direct imaging of the inner 
milli-parsec of a quasar would require angular resolutions on the 
order of micro- or even nano-arcseconds. There is currently no 
instrument capable of such performance. Fortunately, nature has 
provided us with a cosmic magnifying glass in the form of grav- 
itational lensing, which helps to resolve the source, and to reveal 
the spatial structure of its most inner parts. 

Our target is QSO 2237 -f 0305. It was discovered by 
iHuchra et al.l d 1985 1) during the Center for Astrophysics Redshift 
Survey, and is also known as "Huchra's lens" or the "Einstein 
Cross". It is probably the most studied gravitationally lensed 
quasar It consists of a = 1 .695 source gravitationally lensed 

* Based on observations made with ESO Telescopes at the Paranal 
Observatory under program ID 073.B-0243, 074.B-0270, 075.B-0350, 
076.B-0197, 177.B-0615. 



into four images arranged in a cross-like pattern around the nu- 
cleus of a z/ - 0.039 baiTed Sab galaxy. The average projected 
distance of the images from the lens center is 0.9". 

Gravitational lensing is achromatic, in the sense that pho- 
tons are deflected the same way regardless of their wavelength. 
However, lensing is sensitive to the source size, and quasar ac- 
cretion disks are known to vary chromatically from the center to 
the edge. As a consequence, microlensing of the disk by stars 
located in the lensing galaxy will affect the quasar images inde- 
pendently, in ducing chromatic difference s in the spectra of the 
images (Wa mbsganss & Paczvriskilll99ll) . These differences are 
directly linked to the energy profile of the disk, i.e. the scaling 
between the wavelength and the coiTesponding size of the source 
emitting at that wavelength. As a consequence, microlensing- 
induced flux variations can be used to constrain the energy pro- 
file of the quasar accretion disk . 

' Rauch & BlandfordI (1199 Ih were the first to use this tech- 
nique, and for the Einstein Cross, they found that the near- 
ultraviolet emission regions were smaller than expected for ther- 
mal emission from an accretion disk. However, more modern 
work in this system and other gravitational lenses have gener- 
ally reached the opposite conclusion. For instance, ' Poolev et al.l 
(2007) finds that the near-ultraviolet emission regions of ten 
lensed quasars are a factor of 3 - 30 larger than the size pre- 
dicted by simple accretion disk models to produce the observed 
optical flux. In a more quantitative analysis. lMorgan et al.l(l2007l) 
use the microlensing variability observed for nine gravitationally 
lensed quasars to show that the accretion d isk size is consistent 
with the expectation from thin disk theory jShakura & SunvaevI 
Il973l) . However, these sizes are larger, by a factor of ~ 3, than 



2 



A. Eigenbrod et al.: Microlensing variability in the Einstein Cross 



the size needed to produce the observed infrared flux by thermal 
radiation from a thin disk. 

The Einstein Cross is particularly well suited for microlens- 
ing studies, because of the symmetric configuration of the lensed 
images, which results in very short time delays, and because of 
the low redshift of the lensing galaxy, which places the lensed 
images in the bulge of the lensing galaxy, where the probability 
of microlensing by stars is high. Because of this, the Einstein 
Cross has been intensely monitored by di f ferent projects (e.g. , 
Corrigan et al. lll99ll: lOstensen et alJ[T996l: [Alcalde e'tani2002t 
Schm idt et all l2002h . The project having the longest duration 
and the best sampling rate is th e Optical Grav itat ional Lensing 
Exper iment (OGLE) CWozniak et al.l l2000allbt lUdalski et all 
120061) . OGLE has monitored QSO 2237+0305 since 1997, and 
delivers the most complete lightcurves for this system. 

Unfortunately, these monitoring campaigns are usually con- 
ducted in one photometric band, and hence they can be used 
to constrain the source size, but not the energy profile, which 
requires at least two bands. Multi-band monitoring is one ap- 
proach, but even more effective is long-term spectrophotomet- 
ric m onitoring, as described in our first paper (Eigenbrod et al. 
12008 , paper I in the following). In this previous contribution, we 
describe our observations and data analysis, and report signifi- 
cant flux variations in the continuum and broad emission lines 
of the spectra of the four lensed images, indicating that both 
the continuum emitting region and the broad line region are mi- 
crolensed. 

Quasars are known to vary intrinsically, and in order to ex- 
tract the microlensing-induced flux variations, we need to re- 
move the intrinsic fluctuations from the lightcurves. In practice 
we do this by considering the difference between the light curves 
of two images. In the present study we focus on the lightcurves 
of images A and B, because they are the two images undergo- 
ing the strongest flux variations within the time span of our ob- 
servations. In the case of QSO 2237+0305, we can neglect the 
time delays between the images, because th ey are expected to be 
on th e order of one day or less (S chneider et al.ll 19881 : iRix et alj 
119921: IWambsganss & Paczvriskilll994 ). We also know that the 
observed flux F of a lensed quasar image is the product of the 
unlensed flux Fq of the source, the extinction e"^, and the macro 
and micro-magnification fimacm and /v,„,cro 

F — f^niicro f^macro ^ . 

As observed in paper I, image A is affected by long-term (more 
than 5 years) microlensing. Hence the micro-magnification 
fJ^micro has two Components over the time span of our observa- 
tion: one constant long-term magnification fimkro, and one vari- 
able short-term /i;„;^^^(0 

H- micro (0 = 

We shall mention that this separation in long and short-term 
microlensing is purely empirical, and only reflects the impos- 
sibility of retrieving the intrinsic source flux. The difference 
in magnitude Am between image A and B is given by Am - 
-2.5 logipA/ Fb). As observed in paper I, the extinction remains 
constant in time in all four images of the Einstein Cross. The 
macro-magnification given in Table [3]results in a constant mag- 
nitude difference on the order of -0. 1 mag between images A 
and B. As a consequence, the time variability that we observe in 
Am is only due to short-term microlensing 
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Table 1. Journal of the second part of our spectroscopic moni- 
toring of QSO 2237+0305. The JuUan dates are given in HJD- 
2450000. 



ID 


Civil Date 


HJD 


Mask 


Seeing ["] 


Airmass 


32 


10-07-2007 


4292 


1 


0.63 


1.153 








2 


0.55 


1.132 


33 


15-07-2007 


4297 


1 


0.57 


1.158 








2 


0.54 


1.220 


34 


25-07-2007 


4307 


1 


0.79 


1.231 








2 


0.83 


1.161 


35 


03-08-2007 


4316 


1 


0.68 


1.412 








2 


0.68 


1.278 


36 


27-08-2007 


4340 


1 


0.88 


1.133 








2 


0.89 


1.143 


37 


06-09-2007 


4350 


1 


0.67 


1.396 








2 


0.72 


1.252 


38 


20-09-2007 


4364 


1 


0.73 


1.230 








2 


0.83 


1.167 


39 


23-09-2007 


4367 


1 


1.23 


1.530 








2 


1.10 


1.344 


40 


05-10-2007 


4379 


1 


0.59 


1.153 








2 


0.55 


1.132 


41 


10-10-2007 


4384 


1 


0.76 


1.283 








2 


0.64 


1.195 


42 


15-11-2007 


4420 


1 


1.07 


1.189 




16-11-2007 


4421 


2 


0.80 


1.148 


43 


01-12-2007 


4436 


1 


1.00 


1.502 



2 0.82 1.318 



Table 2. Wavelength intervals and central wavelength (in A 
measured in the source frame) of the six photometric bands. 



Band 




A, 


1 


[1500, 1750] 


1625 


2 


[1750, 2000] 


1875 


3 


[2000, 2250] 


2125 


4 


[2250,2500] 


2375 


5 


[2500, 2750] 


2625 


6 


[2750, 3000] 


2875 



where mo is a constant. Microlensing depends on the source 
size. Smaller sources are more strongly affected by microlens- 
ing than larger sources. We also know that bluer photons are 
emitted closer to the center of a quasar than redder ones. 
Hence we expect stronger variatio ns of Am at bluer wavelengths 
(Wambsganss & Paczvriskilll99Ih . 

In this paper, we use this chromatic behavior of microlens- 
ing to compare the observed variability of Am at different wave- 
lengths with numerical simula tions. Using B ayesian analysis 
similar to Kochanek (2004b and lAnguita et al.l (|2008), we derive 
the probability distributions for the size of the source emitting at 
a given wavelength. This eventually defines the energy profile of 
the quasar accretion disk. 

In the following, we consider a flat cosmology with 
(Q„, Qa) = (0.3, 0.7) and Hq = 70 km s"' Mpc"'. 

2. Observations 

We use two different data sets for this study. The first one 
is the well sampled V-band lightcurves of the OGLlJ] project 

' http://www.astrouw.edu.pl/~ogle/ogle3/huchra.html 
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Table 3. Macro -lensing parameters for images A and B from 
lKochanekl(l2004l) 



Image 


K 


J 


f-^inacro 


A 


0.394 


0.395 


4.735 


B 


0.375 


0.390 


4.192 



j Udalsldet al1l20Q6l) . from which we select the data between 
Julian days 2453126 (April 30, 2004) and 2454439 (December 
4, 2007), and which comprises 352 data points. We bin the data 
points separate by less than one day, and are left with 181 points 
for which we recompute the errors as described in section |4] 

The second data set is our deep spectrophotometric mon- 
itoring obtained with the FORSl spectrograph mounted on 
the Very Large Telescope (VLT) of the European Southern 
Observatory (ESO). In paper I, we presented the spectra of all 
four quasar images on 3 1 different epochs between October 2004 
and December 2006. We have now completed one more year of 
monitoring, and a total of 43 different epochs spanning more 
than three years until December 2007 are now available. The 
journal of the new observations is given in Table [1] We remind 
that we use two masks for the spectroscopic observations. The 
mask 1 has a slit centered on quasar images A and D, while the 
mask 2 has a slit centered on B and C. In the present study, we 
are interested in the difference magnitude Am between images 
A and B. Thus, we will only consider the epochs where both 
masks 1 and 2 (see paper I), and hence both quasar images have 
been observed on the same night. We therefore drop the epochs 
with ID 7, 8, 10, and 30, which reduces the number of usable 
epochs to 39. The reduction and further processing of the new 
data are done following the same procedure as described in pa- 
per I; we spatially deconvolve our spectra to remove the lens- 
ing galaxy, and we decompose the spectra into a sum of broad 
emission lines, continuum and iron pseudo-continuum. We fit 
the continuum with a power-law 

where Aq - 2000 A in the source frame, and where ay is the 
commonly used canonical power index. 

Based on this method, we could identify strong variations in 
the exponent and intensity /o in the spectra of quasar images 
A and B. This indicates that the amplitude of the variations are 
significantly higher in the blue than in the red part of the spectra. 
In order to study in more detail these chromatic fluctuations of 
the continuum, we define six wavelength ranges, that we will use 
as photometric bands. The bands all have the same width and 
cover the whole wavelength range between 1500 and 3000 A 
measured in the source frame, see Table |2] For each of these 
bands we compute the Am lightcurve by integrating the contin- 
uum power-law in the corresponding wavelength range. We give 
the corresponding AB magnitudes of the four quasar images in 
Table |6] By subtracting the magnitude of image B from that of 
image A, we obtain six 39-epoch lightcurves, as displayed in 
Fig. [T] We immediately see that the lightcurves for the redder 
bands (e.g. band #6) vary less that the bluer ones (e.g. band #1). 

3. Microlensing Simulations 

In order to simulate the microlensing effects, we have to choose 
a macro-model for the lensing galaxy. Many models have been 
proposed for the Einstein Cross (see for instance the summary in 



variability in the Einstein Cross 3 



1 




\ I I I I I 

3500 4000 
HJD-2450000 



Fig. 1. The difference lightcurve Am between quasar images 
A and B measured in different wavelength bands. For clarity 
we shift the curves along the vertical direction. The binned 
OGL E observations with 181 data points measured in the V- 
band dUdalski et al.l2006h . From our 39-epoch spectrophotomet- 
ric monitoring, we compute the difference lightcurves of the con- 
tinuum integrated in the six bands defined in Table|2l 



Table 2 of'Wyithe et al."2002'). We choose the macro-model pro- 
posed by Kochanek (2004). The surface density k, the shear y, 
and the corresponding macro-magnifications jimacro are given in 
Table [3] The convergence k is usually separated into a compact 
matter distribution a-* and a smooth matter distribution Kc com- 
posed of gas and dark matter. However, as the quasar images are 
located within the bulge of the lensing galaxy, we do not expect 
a strong contribution from the smooth matter tKochanek et al.l 
2007), so we neglect Kc, which amounts to assuming the entire 
surface mass density is in stars (i.e. k - /<•*). 

The next step is to define the mass distribution of the mi- 
crolenses. Earlier investigations of microlensing have led to 
the conclusion that the magnification probability distribution 
is not very sensitive to the mass function of the microlenses 
(e.g. Wambsganss 1992; Lewis & Irwin 1995; Wvithe & TurneB 
l200ll ; ICongdon et alj|2007l) . but that it depends on the mean mi- 
crolens mass (M), which sets the parameter scales, such as the 
Einstein radius. The determination of (M) is not trivial because 
of the degeneracies existing between (M), the effective trans- 
verse velocity, and the source size. Several studies succeeded in 
getti ng an estimate of (M), and currently favor (M) ^ 0.1 Mq 
(e.g. IWambspanss et alJ[T990bt iLewis & Irwinlll996l) . We will 
assume this value for all the microlenses distributed in the lens 
plane. This value is also consistent with the mean mass found by 
the MACHO collaboration who studied microlens ing from ob- 
jects located to wards the bulge of the Milky Way dAlcock et al.l 
Il997). See also IWvithe e"taiT ( l2000bl) . who give a summary of 
the determinations of microlens masses in other galaxies. 

The masses of the microlenses sets the Einstein radius r^. 
For QSO 2237-1-0305, the Einstein radius projects in the source 
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plane as 



rE = D, = 5.77 X 10>^ (M/O.lMo)'/^ cm. 

Having defined k, y and the mass of the microlenses, we can 
go one step further in our simulations, and start generating large 
magnification patterns for the two quasar images A and B. We do 
this u sing the inverse ray- shooting method (Wambsganss 1990; 
IWambsganss et al]|1990al) . For each image, we shoot approxi- 
mately lO" rays from the observer through the lens plane back 
to the source plane, where the rays are collected in a 10, 000 by 
10, 000 pixels array, corresponding to 100 by 100 Einstein radii. 

In order to study a sample of different source sizes, the mag- 
nification patterns are convolved with a set of source profiles. 
Microlensing-induced flux variations are relatively insensitive 
to the source profile, b ut they depend on its characteristic size. 
iMortonson et al.l (|2005) showed that the half-light radius is the 
most important quantity for controlling the shape of microlens- 
ing light curves, whereas the source profile is less important. For 
simplicity we choose a Gaussian profile for the surface bright- 
ness, and we define the size of the source as the Full Width at 
Half Maximum (FWHM), which is twice as large as the half- 
light radius, and 2.35 greater than the Gaussian width cr. 

The relative size of source with respect to the Einstein radius 
controls the smoothing of the magnification pattern. Very large 
source sizes (i.e. larger than 4 r^) produce magnification patterns 
that are so strongly smoothed, that they are unable to account 
for the amplitude (i.e. higher than 1 mag) of the flux variations 
seen in the O GLE hghtcurves. Hence, as was sho wn in previous 
studies (e.g. Wvithe et al.ll2000atlYoneharal2001 ), we can safely 
rule out continuum source sizes larger than 4 r^. On the other 
hand, very small source sizes (i.e. smaller than 0.01 r^) give 
very moderately smoothed magnification patterns, which lead to 
sharp magnification events that are not observed in the OGLE 
lightcurves. The lower bound found in that way for the source 
size is 0.01 (e.g. Ko chanek 200 4; Anguita et al,.2008.) . 

In the present study, we vary the FWHM of the Gaussian 
profile from 1 to 400 pixels, which corresponds to 0.01 and 
4.00 rE, respectively. We consider 45 different source sizes be- 
tween these two extreme values. They are 1, 2, 4, 6, 8, 10, 12, 
14, 16, 18, 20, 22, 24, 26, 28, 30, 35, 40, 45, 50, 55, 60, 65, 70, 
75, 80, 85, 90, 95, 100, 120, 140, 160, 180, 200, 220, 240, 260, 
280, 300, 320, 340, 360, 380, and 400 pixels. 

4. Microlensing simulations fitting the OGLE data 

O ur simulations are conducted following the method described 
in lKochanekl (|2004|) . This technique is based on a Bayesian anal- 
ysis, that determines the probability distributions for the physical 
parameters of interest by comparing trial lightcurves with the ob- 
served data. Because we want to infer probability distributions, 
we need to simulate a large number of these trial lightcurves 
in order to obtain a statistically significant sample. Thus, for 
each of our 45 source sizes, we simulate 10,000 light curves 
for both quasar images A and B by tracing source trajectories 
across the corresponding magnification patterns. We extract the 
pixel counts along the positions of the trajectory using bi-linear 
interpolation and we convert them into magnitudes. We subtract 
the simulated lightcurve of image B from that of image A, and 
obtain a total of 4.5 x 10^ diff'erence Hghtcurves (10,000 for 45 
different source sizes). Eventually, we want to compare this large 
library of simulated difference lightcurves to the Am observed in 



file OGLE data, but before we do that, we specify some charac- 
teristics of the simulated source trajectories. 

The simulated lightcurves are obtained from source trajecto- 
ries across the magnification patterns, and these trajectories are 
characterized by parameters of two kinds: physical and trajec- 
tory parameters. The physical parameters are the local magni- 
fication tensor k and y, the mass (M) of the microlenses, the 
size Rs of the source, the effective transverse velocity V of 
the source, and the magnitude offset mo between the two im- 
ages due to a combination of the macro-magnification, long- 
term microlensing, and differential extinction between the im- 
ages. Because the scales of the magnification patterns are de- 
fined in terms of Einstein radii, our computational variables 
for the source size and velocity are in fact the scaled source 
size Rs - Rs/{M/Q.l M©)'^^ and the scaled transverse veloc- 
ity V - V/{M/0.l M^y^^. Foflowing these definitions, Rs is the 
relative size of the source with respect to the Einstein radius, 
and V is the velocity with which the source is moving across the 
magnification pattern. 

Each source trajectory is also defined by the trajectory pa- 
rameters (0, xo,a,xo,b), where is the direction angle, xo,a are 
the coordinates of the starting point in the magnification pattern 
of image A, and xq.b in the pattern of image B. The trajectory is 
constrained to have identical directions 6 and velocities V across 
the patterns of images A and B. The direction is set to be the 
same in both patterns because the shear di rection betwe en im- 
ages A and B is approximately the same (IWitt&MadlI994 . 
and because we assume the motion of the source to be primarily 
due to the bulk motion of the lensing galaxy rather than to th e in- 
dividual motions of the stars. Kundic & Wambsganss! (ll993l) and 
Wambsganss & Kundic ( 1995) show that the velocity dispersion 
of the stars can statistically be interpreted as a bulk velocity ar- 
tificially increased by an efficiency factor a ^ 1.3. 

For each particular choice of the parameters, that we write 
p = {k, y, (M), Rs, V, mo, 0, xq^a, xo,b), we get one simulated dif- 
ference lightcurve Ain'j^ip), that can be compared to the observed 
data D - Am^ by measuring the goodness of fit with ax^ statistic 

= ) 

where cTk are the uncertainties of the OGLE data, and nohs = 181 
is the number of the binned OGLE data points. We determine 
the (Til from the photometric errors ctogle given by the OGLE 
project, from the standard deviation between the binned points 
CTbin, and from the systematic error of the OGLE data ctj,,,, with 

~ '^OGLE '^hin '^^vs- We estimate the systematic error 
cr,,,, by carrying out a polynomial fit of the OGLE difference 
lightcurve Anik, using a polynomial of high enough order that 
the residuals are uncorrected, i.e. such that the auto-correlation 
function of the residuals reduces to less than 0.5 within a sepa- 
ration of a few data points. We find that a polynomial on order 7 
is sufficient, and ctj,., = 0.03 mag. 

By varying the parameters we construct our library of 4.5 x 
10^ lightcurves. This produces an ensemble of models, which 
given the data and using Bayesian analysis (see Section |5]l can 
infer the probability distributions for the parameters. The size of 
the magnification patterns is approximately one hundred times 
greater than the time scale of the observations multiplied by the 
effective velocity. Hence the available parameter space is huge, 
and it is relatively easy to find good fits to the data. 

In our simulations we fix the values of k, y, and (M), 
and are left with the following set of variable parameters p = 
{Rs, V, mo, 6, xo,A, Xo,b). The physical parameters that we want to 
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data D - Amt, is 
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Fig. 2. Examples of six simulated lightcurves fitting the observed 
OGLE difference lightcurve Am between quasar images A and 
B. 



investigate are the source size and the effective transverse veloc- 
ity. The other parameters, i.e. (mo, 0, xo,a,xo,b), have no direct 
physical implications, and can be chosen arbitrarily. They are 
therefore called nuisance parameters. The influence of these pa- 
rameters on the inferred probability distributions of Rs and V 
vanishes when using Bayesian analysis and a sufficiently large 
library of simulated lightcurves, as described in Section |5] 

A first guess for the nuisance parameters (mo, 0, X().a,xo,b) 
is chosen randomly following a uniform distribution. In the 
Bayesian analysis, lightcurves with low Xqcle '^^lues have rela- 
tively higher weights, and contribute more to the final probabil- 
ity distributions. An effective and fast method consists in search- 
ing for trajectories with ^ow x^qie values. To do this and hence 
minimize the nec essary comput i ng tim e, we follow the proce- 
dure described in lAnguita et al.l (l2008l) . and optimize the seven 
parameters (V, mo, 9, xo,a, xo,b) with a;^^^ -based minimization al- 
gorithm using a Levenberg-Marquardt least squares routine in 
IDI0. We verify that this minimization technique still samples 
uniformly the whole magnification pattern. Finally we obtain a 
trajectory library containing 4.5 x 10^ trajectories (10,000 tra- 
jectories for 45 different source sizes) fitting the OGLE data. All 
the source trajectories have a reduced Xqqh; /nd,,/ < 10, where 
ndof = 345 is the number of degrees of freedom, i.e. the number 
of observations minus the number of fitted parameters (here 7). 
Six examples of the simulated lightcurves are given in Fig.|2] 



5. Bayesian analysis 

In the previous section, we have described how we build 
our large library of trial lightcurves. We will now use these 
lightcurves to determine the probability distributions for the pa- 
rameters Rs and V based on Bayesian analysis as described by 
iKoch anek ( 2004). Following Bayes' theorem the probability of 
the parameters /5 = (/c, y, (M),^s, V, mo,0, Xo,a,xo,b), given the 



Pip\D) = 



P(D\p)P{p) L(D\p) 



P(D) 



N 



Pip) 



where P{p) is the prior, L{D\p) is the likelihood and is a nor- 
malization constant. 

We do not have any prior information either on the trajectory 
parameters 0, xo,/\,Xo,b or on the magnitude offset mo. Therefore 
we choose these parameters from a random and uniform distri- 
bution. This means that the corresponding priors P(mo), P{9), 
f (xo,a), and P(xo,b) are constant. We also assume that the pa- 
rameters are independent, i.e. that 

P{p) = P{k) P(y) P{{M)) P(Rs) P(V) P(mo) P{0) P{xq^) P(xo,b) ■ 

We define the relative likelihoods of the parameters p based on 
the^^ statistics. Usually, this is d one following the standard ap- 
proach for ensemble analysis (e.g.. lSambridgd 19991) . which uses 
the maximum likelihood estimator 



L(D\p)^exp{~x\p) 



However, as noted by Kochanek (2004), using this standard esti- 
mator for the OGLE data works poorly because we are compar- 
ing the probabilities of completely different light curves rather 
than models related to each other by continuous changes of pa- 
rameters. T o circumvent this , we use the likelihood estimator 
proposed bv lKochaneS (|2004 



ndof-2 xHp) 



UD\p) = r 



where F is an incomplete Gamma function. Finally, all Bayesian 
parameter estimates have to be normalized by the requirement 
that the total probability is unity, i.e. J P(p\D) dp = 1 and hence 



j L(D\p)P(p)dp. 

In practice we sum the probabilities for our random sampling of 
trajectories, which is equivalent to using Monte Carlo integra- 
tion methods to compute the integral over the ensemble of all 
possible trajectories. The sum over the random trajectories will 
converge to the true integral provided we consider a sufficiently 
large number of trajectories. The probability of a trajectory j 
(defined by the set of parameters pj) given the data D is 



P(pj\D) 



L{D\pj)P{Rs,j)P(Vj) 
Y:j^yL{D\pj)P{Rsj)P{Vj) 



(1) 



http://cow.physics.wisc.edu/~craigm/idl/ 



where n is the total number of trajectories in the library, and 
P{Rs,j) and PiVj) are the priors on the source size and velocity, 
respectively. The priors on the other parameters are the same for 
all trajectories, and cancel out in the fraction above. The final 
probability distribution of the parameters are obtained by sum- 
ming the probabilities P{pj\D) of all the trajectories j of the li- 
brary. 

The prior P{Rs,j) is required because we have not uniformly 
sampled the parameter space of the source size. We have to cor- 
rect for this. The prior P{Rs.j) is related to the density of tra- 
jectories as a function of R^j given our sampling of the source 
size in 45 different bins. Each bin contains 10,000 trajectories. 
Therefore, small bins have higher densities, and the correspond- 
ing PiRsj) should be proportionally lower to compensate this 
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IQis 1015 IQi^ 1000 10' 10^ 13 3 4 

R, ({M/0.1 Mg)'^^ cm) V ((M/0.1 U^Y'^ km/s) R./V (10'^ cm/(km/s)) 



Fig. 3. Left : The probability distributions for the FWHM ^.5 of the source, based on the OGLE data. We consider two cases, one 
with and another without the prior on the velocity. Center : Same as left panel, but for the effective transverse velocity V. Right : 
The probability distributions for the ratio Rs/V- 



effect, i.e. if the size of bin b is lb, and if falls in this bin, 
then P(Rsj) oc Thus each bin has the same density of trajec- 
tories. 

As already mentioned, there is a strong, essentially linear, 
correlation between R^ and V, i.e. R^ <x W where x ^ I 
(Kochanek 2004). Furthermore, Rg and V are both defined as 
functions of (M), which imply degeneracies between (M), V, 
and Rs. However, since Rs = ^^(M/O.lMo)'/^ and (M) cc 
(y/y)2, the physical size of the source Rs oc y y * ' ^ y de- 
pends essentially on our estimate of the physical velocity V, and 
avoids the degeneracies between (M), V, and Rs- 

Estimates of V can be obtained from the obser vations of the 
motion of other galaxies or galaxy clusters (e.g. ' Benson et al.l 
l2003h . which show that peculiar velocities are typically not 
higher than 1, 500 km/s. In order to give a lower weight to tra- 
jectories having much higher velocities, we have to introduce 
a prior on the velocity. Previous studies have considere d vari- 
ous priors (e.g. lKochanei3l2004l: iGil-Merino et al. I l2005h . some 
are more restrictive than others, but most of them have in com- 
mon that they favor transverse velocities of ~ 600 km/s in the 
lens plane, which correspond to projected transverse veloci- 
l .ies of ~ 6,000 km /s in the source plane of QSO 2237H-0305 
jKavser et al.lll986 l). 

The effective transverse velocity is the result of the rel- 
ative motion betwe en the source, the lens and the observer 
dKavser et al.l 1 19861) . and is enhanced by a contribution from 
the velocity dispersion of the stars in the lensing galax y 
(|Kundic & Wambsgansi [19931 IWambsganss & Kundid Il995l) . 
The problem is that we do not known the peculiar velocity of 
either the source, or the lens. The best we can do is to consider 
probability distributions for these unknown velocities, based on 
what we know from the peculiar motion of other galaxi es (e.g. 
iMould et al.lll993l:lBenson et al.ll2003h . kochan ek (2004) shows 
that a good approximation for the probability distribution of the 
effective transverse velocity in the source-plane is 



Table 4. Results from the microlensing simulations and 
Bayesian analysis applied to the OGLE data. 



y 

PiV) = /„ 



yy 



cr- 



exp 



y2 + y2 

2cr2 



(2) 



where Iq is a modified Bessel function, V ^ 2500 km/s is the 
mean velocity, and cr ^ 3500 km/s is the total root mean square 
velocity obtained by considering all the different velocities con- 
tributing to the transverse effective velocity. In our analysis we 



with velocity prior 

= (9.2+«;») X 10'^ <M/0.1Mo>'/2 cm = (0.16+;;;|5) 

P = (3.9+]-«) X 10^ <M/0.1Mo)'/2 km/s 

RJV = (2.1 ± 0.5) X 10'- cm/(km/s) 
without velocity prior 

R, = (4.0:^;^) X 10"^ (M/0.1Mo>"2 cm = (0.69:^;^°) re 

V = (l.8+2-^) X 10* <M/0.1Mq>i/2 km/s 

R,JV = (2.1 ± 0.6) X 10'2 cm/(km/s) 



assume (M) = 0.1 Mq, and we define the prior on the scaled 
velocity V as P(V) = P(V). 

Using Bayesian analysis, and considering both cases with 
and without velocity prior, we obtain the probability distribu- 
tions of the source size ^5, the effective transverse velocity V, 
and the ratio Rs/V plotted in Fig. |3] The results are given in 
Table g] 



6. Energy profile of the quasar accretion disk 

Accretion disk models for quasars generally make the as- 
sumption that the disk is optically thick and geometrically 
thin; both assumptions are required to cause efficient conver- 
sion of the gr avitational potential energy of accreting matter 
into radiation (Pringle & Rees" 19721: IShakura & Sunvaevll 19731: 
Novikov & Thorne 1973). Alternative models in which mass is 
lost in a wind dKuncic & Bi cknell 2007), the disk accretes at 
greater than the Eddington limit and thus becomes geometrically 
thick, in which case adv ection or convection can dominate the 
radial heat transport (e.g.'Abram owicz et al.lll98 "8). or magnetic 
connections between the black ho le and disk modify the rate of 
energy dissipation (lAgol & Krolikll2000l) . 

The peak of the continuum emission in a typical quasar spec- 
tral energy distribution is in the "Big Blue Bump", and can be 
fitted well with a spectrum which scales as Fy oc vi/3g'"'/*r^ 
with a characteristic temperature T -1 x K (Krolik 1998lJ). 
This is predicted by the standard thick accretion disk model in 
which r oc /;-3/4( 1 _ {R.^ IE) 1/2)1/4^ the functional form predicted 
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Fig. 4. Top : The source FWHM ratio Ri/Rref - Ri/Rref as a function of the wavelength ratio Ai/Amf obtained from the trajectories 
in the library having R, = 0.1 - 10 pixels (left) and Rs - 1 - 100 pixels (right). The reference wavelength is Ar^f = 2125 A 
measured in the source frame. The error bars give the Icr uncertainties. The solid line indicates the linear regression across the 
points. The dashed line shows the relation expected for the standard R oc A'^^^ scaling. Bottom : The corresponding probability 
distributions of the scaled source sizes R, obtained in the six bands / = 1,2, 6, and using the same color code as in Fig.[T] 



by non-relativistic thin accretion disks with an inner edge Rj„ 
at which the viscous stress disappears (relativity modifies this 
somewhat - lKrolik|[T998ai) . At large radius T oc R^^l^, and the 
model predicts that the radiative flux F (erg cm"^ s ') at the 
surface of the disk at a radius R is proportional to the energy 
production rate at that radius, resulting in F oc R^^. Since we are 
assuming a blackbody spectrum, F oc T"*, and because the wave- 
length A of the peak of the black body spectrum is proportional 
to r^', most of the radiation at a given radius R comes from near 
a radius 

R oc T-''' oc ^4/3 . 



This defines the energy profile of the quasar, away from the in- 
ne r edge of the d isk, which is identical to the form obtained 
bv lKochanek et al.l (i2007l) from the standard thin accretion disk 
model of lShakura & Sunvaevl ( ll973l) . 

However, given that we are trying to test the accretion disk 
model, we allow the temperature to scale as a power-law with 
radius with arbitrary slope T oc R^^ (R auch & Blandford 199 1|) ; 
this functional form is general enough to cover a variety of dif- 
ferent alternatives to the standard accretion disk model. 

There exist numerous models for quasar accretion disks in 
the literature. We can not consider all of them here, and we se- 
lect only a few presenting an interesting variety. One particular 
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model is motivated by the fact that the theoretical optical/UV 
continuum of a standard thin disk is Fy oc v^/^, which is in- 
consistent with Fy oc V ' observed in many quasars. Hence, 
iGaskelll (2008) suggests that the observed quasar spectra can 
be reproduced by accretion disks with a temperature gradient 
of r oc instead of T oc R^^^'^. This implies 



Roc A 



1.75 



Another model is proposed by lAgol & KrolikI ( |2000|) . and de- 
scribes a disk which is powered by the spin of the central black 
hole. This changes the temperature profile to T oc r"^^*^, and pre- 
dicts 

R oc //^ . 

The three models mentioned above all predict a power-law rela- 
tion between the wavelength and the corresponding source size 

RozA^ 

with different values for the index ^. For two different radii R 
and Rref emitting at the wavelengths A and A^ef respectively, we 
have 



R 

Rref 



A 

Aref 



log 



R 

Rref 



^ log 



Aref 



To distinguish between these different accretion disk models, 
we determine in the next section, which power-law is best com- 
patible with our spectroscopic observations, and which of these 
models can be ruled out for QSO 2237H-0305. 



7. Interpretation of the spectroscopic data 

7.1. Method 

We will now compare the theoretical energy profile with our ob- 
servations. In the procedure described in the previous sections, 
we have built a library of simulated lightcurves by tracing a large 
number of source trajectories through the magnification patterns 
for various source sizes in order to reproduce the OGLE V- 
band photometry Am of QSO 2237-1-0305. We use this ensemble 
of models as an input to extend our analysis to our spectroscopic 
observations. We consider the integrated lightcurves in the six 
color bands defined in Table|2]as shown in Fig.[I] For each band 
i, we have a lightcurve Am(/},, f) with 39 epochs. The wavelength 
Ai is the central wavelength of band / as defined in Table |2] We 
then determine which scaled source size R, best reproduces the 
observed lightcurve in band /. 

Each source trajectory /' in the library is defined by a set of 
the parameters pj = (Kj,yj,{M}j,R,j,Vj,moj,0j,XQ^A,j,TiQ,B,j)- 
For every trajectory j, we keep all these parameters fixed, ex- 
cept the source size Rsj, that we modify, i.e. we trace the same 
source trajectory in all the convolved magnification patterns. By 
doing this, we get the simulated difference lightcurve Am'j(t) cor- 
responding to the 45 different source sizes defined in Section [3] 
We then interpolate the magnitude values of these extracted light 
curves to obtain continuous values for the source size. For ev- 
ery trajectory j, this defines a surface Am'j{Rs, t) in the three- 
dimensional space defined by Rs, t and Am. 

We use these Am'j(Rs, t) to interpret our spectroscopic data 

by determining which scaled source size Ri best reproduces the 
observed lighcurve Am(/l,, f) in band The size Rj is the value 
that minimizes the x]j between the spectrophotometric observa- 
tions AmiAi, t) and the simulated difference lightcurve Am'^iRj, t). 



Gaskell (2007) 



"Shakura & Sunyaev (1973) 



0IT& Krolik (2000) 



R. ({M/0.1 Mg>'/2 cm) 

Fig. 6. The power-law index ^ of the energy profile oc /l^ as a 
function of the scaled FWHM Rs of the source. The points are 
obtained from the subsamples of trajectories with given R^ fit- 
ting the OGLE data and including the prior on the velocity. The 
horizontal solid red lines indicate the index ^ and its Icr devia- 
tion obtained by considering the whole library of trajectories and 
including the velocity prior. The dashed lines show the expected 
indices for the three indicated accretion-disk models. 



After this second fitting procedure, we update the relative likeli- 
hoods of the lightcurves 



L{D\pj,i)^Y 



X exp 



2 \ 



where / refers to the considered band. In the fitting of the 
spectroscopic observations Am{Ai, t), we are comparing similar 
light curves, where only Ri is changing. This is a very differ- 
ent situation than when we were fitting the OGLE difference 
lightcurve and comparing simulations with very different param- 
eters. Hence, we can use here the standard maximum likelihood 
estimator exp(-xjj/2) to compute the likelihoods L{D\pj, i). We 
run this procedure for all the simulated source trajectories j in 
the library, and use these updated likelihoods to compute the fi- 
nal probability distribution for the parameters of interest, as de- 
scribed in section |5] 

7.2. Results 

For each band /, we compute the probability distribution of the 
source-size ratio Rj/Rref following the Bayesian analysis de- 
scribed in Section |5] using equation ([T]i. The reference radius 
Rref at the reference wavelength Aref can be chosen arbitrar- 
ily, and we simply define the middle band #3 as our reference. 
The scaling between R, and ^ 5 is independent of the considered 
wavelength and Ri/Rref = Rj/Rref- Thus, the ratio Rj/Rref is 
not expected to depend on the assumed microlens mass (M). We 
plot the ratio Rj/Rref against the corresponding wavelength ratio 
Aj/Aref, and determine the slope of the energy profile by fitting a 
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Fig. 5. Left : The final energy profile obtained from the whole ensemble of source trajectories, including the prior on the velocity. 
The reference wavelength is Are/ = 2125 A measured in the source frame. The error bars give the Icr uncertainties. The solid line 
indicates the linear regression across the points. The dashed line shows the relation expected for the standard R oc A'^^^ scaling. Right 
: The corresponding probability distributions of the Ri values obtained in the six bands / = 1,2, 6, and using the same color code 
as in Fig.[T] 



Table 5. The relative size Ri/Riej obtained for the six photomet- 
ric bands, including the velocity prior. 



Band 


Ac [A] 


R.IRref 


1 


1625 


0.73 ±0.14 


2 


1875 


0.87 ±0.11 


3 


2125 


1.00 ±0.09 


4 


2375 


1.13 ±0.09 


5 


2625 


1.29 ±0.13 


6 


2875 


1.46 ± 0.18 



power law 

RocA^. 

We do this for the entire library of source trajectories, but also 
individually for every subsample of the library with the same 
initial source size R^ used to fit the OGLE data. The resulting 
index f is plotted as a function of Rs in Fig. |6] We give two 
examples of energy profiles in Fig.|4]derived from subsamples of 
trajectories: one having Rs -QA fe - 10 pixels and another with 
Rs - I te - 100 pixels. The final energy profile obtained from 
the whole sample of source trajectories is shown in Fig.|5] given 
in Table |5] and yields ( - 1.2 + 0.3 when adding the velocity 
prior Without this prior we get a flatter energy profile with ^ - 
1.1+0.3. 

8. Discussion 

From the results of our microlensing simulations, we can ex- 
tract several interesting facts. First, the source FWHM Rs - 
(0.16;^Q [g) we derive from the fitting of the OGLE data (in- 
cluding the prior on the velo city) is wel l compatible with the up- 
per limit of 0.98 given b v I Yoneharal (|200 li) . Our result is also 



in goo d agreement with the FWHM derived by lAnguita et al.l 
( I2008h . Rs = (0.06 + 0.01) r^, and bv lKochanekl (|2004|) . = 
(0.20^Q [2) rE- The differences between these estimates of Rs 
can have two origins. First, as mentioned by iKochanekl (|2004|) . 
the probability distributions obtained for the parameters depend 
on the choice of the period of observations. If the considered 
period is very active in terms of microlensing, the simulations 
will favor high transverse velocities and/or small source sizes. 
Second, the choice of the velocity prior has a strong effect on 
the derived source size. As we observe in Fig. [3] there is a 
strong correlation between the scaled source size Rs and the 
scaled effective transverse velocity V. Independently of the ve- 
locity prior considered, we find that the source size is directly 
proportional to the transverse velocity, with the relation RsjV - 
(2.1 + 0.6) X 10^^ cm/(km/s). This is in excellent agreement with 
IKochanekl (12004 . who also finds RslV = 2.1 x 10'^ cm/(km/s), 
where Rs is expressed as the FWHM. This defines the time 
scale of the observed microlensing-induced fluctuations, which 
is given by the half-light radius divided by the effective trans- 
verse velocity, i.e. 0.5 RslV - 4.0 +1.0 months, and which is 
independent of microlens masses or velocity priors. 

The strong correlation between ^5 and V (lKochanekll2004l) 
implies that the choice of the velocity prior and the selection of 
the observation period can both bias our estimate of the source 
size. The microlensing activity during our observations was not 
as important a s du ring the periods studied_by Kochanek ( 2004 ) 
and lAnguita et alj fcOOSi) . Furthermore, lAnguita et alj (120081) 
conside red a shorter observatio n period, and used the velocity 
prior of iGil-Merino et a (120051) . which is more constraining to- 
wards high velocities (hence favoring smaller source sizes) than 
the one we considered in the present study. This explains, why 
our derived Rs is slightly larger One possible solution to mini- 
mize the bias induced by the selection of the observation period 
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is to consider the longest possible lightcurves, but this has the ev- 
ident drawback of dramatically increasing the computing effort, 
making the problem rapidly intractable. However, as discussed 
in the following, the estimate of the source size has only a lim- 
ited effect on our determination of the energy profile. 

From our library of simulated source trajectories fitting the 
OGLE data, we extend our analysis to our spectroscopic data, 
and determine the energy profile R cc of the quasar accre- 
tion disk. When we consider the whole library and include the 
prior on the velocity, we get = 1.2 ± 0.3, which is in good 
agreement with the model of lAgol & KrolikI ( |2000|) and with 
the ^ = 4/3 index expected from the standard accretion disk 
m odel. It is l ess co mpatible (at Icf) with the accretion disk model 
of iGaskelll (l2008l) . Our findings are also in good agreement 

with lAnguita et all (|2008|) . who obtain (6400 A/4700 A)^ = 

1.45;^Q 2g, which yields a power-law index ( = l-2^0 6' Without 
the velocity prior, we get a lower value ( - 1.1 + 0.3, which is 
still compatible with ^ = 8/7 and 4/3, but less with the model 
of Gaskell ( 2008i) . This is expected because, without the prior, 
large sources gain relatively more importance in the Bayesian 
analysis, and favor flatter energy profile, as shown in Fig. |6] 
However, as discussed earlier, large source sizes imply unlikely 
high transverse velocities, and hence a much more reliable result 
is obtained when including the velocity prior 

We observe that our result for the R cc scaling is not very 
sensitive on the velocity prior with respect to other parameters. 
For instance, we know that the source size ^ 5 is extremely sensi- 
tive to such a prior, and can vary over orders of magnitudes if the 
considered transverse velocity does. On the contrary, the index 
^ varies only by a few percent, even if the considered velocities 
are modified by orders of magnitude. Hence, and independently 
of the determination of V and ^5, the chromatic variations be- 
tween two images of a lensed quasar are extremely efficient in 
constraining the relative sizes of different region s of the accre- 
tion disk. This is further confirmed by the fact that lAnguita et al.l 
(I2OO 8) obtain a value for ^ which is in good agreement with ours, 
even though they derive a slightly smaller source size. 

The influence of the initial value of the parameter ^ 5 (chosen 
to fit the OGLE data) on the resulting energy profile of the accre- 
tion disk is obvious in Fig.|6] Sources having sizes in the range 
0.06 < Rs < 0.2 r£ are in good agreement with the f = 8/7 
and 4/3 scalings, while larger sources give flatter energy pro- 
files, and smaller sources steeper profiles. The same behavior is 
observed, when we omit the velocity prior 

9. Conclusion 

We present the continuation of our spectrophotometric monitor- 
ing of QSO 2237+0305 conducted at the Very Large Telescope 
of the European Southern Observatory, which extends over more 
than three years from October 2004 to December 2007. Our pro- 
gram provides the spectra of the four lensed images at 43 epochs. 
Analysis of these data show that the continuum and the broad 
line region of the background quasar are microlensed, and that 
images A and B are particularly affected. In this paper, we fo- 
cus on the observed variations in the continuum of the spectra of 
these two images, and use them to constrain the energy profile 
of the quasar. 

We build microlensing magnific ation pattern s with 
the in verse ray -shooting method ( iWambsganssI Il990t 
[Wambs ganss et al. 1990a), and convolve them with differ- 
ent source sizes. We trace source trajectories through these 
convolved patterns and fit the observed OGLE difference 



lightcurve between image A and image B. Combining these 
simulations with Bayesian analysis following the method 
proposed by Kochanek (2004), we infer probability distributions 
for the effective transverse velocity V, and the source size R, of 
the quasar Our results are compatible with previous studies. 

Besides the OGLE broad-band photometry, we analyze the 
data of our spectroscopic monitoring and derived the corre- 
sponding difference lightcurves in six different photometric 
bands. Each band is 250 A wide and together the bands cover 
the wavelength range between 1500 and 3000 A (measured in the 
rest frame of the source). We compute the difference lightcurves 
between images A and B in these six bands, and observe that, 
as expected from the microlensing of an accretion disk, bands at 
bluer wavelengths exhibit stronger flux variations than bands at 
redder wavelengths. 

Based on these chromatic variations and using the ensem- 
ble of the microlensing simulations that are fitting the OGLE 
data, we determine the relative sizes of the regions of the accre- 
tion disk emitting in the six photometric bands, i.e. we derive 
the energy profile of the accretion disk. We find that this pro- 
file follows a power-law R o: A^ with ^ = 1.2 + 0.3, which is 
in good agreement with the standard thin accretion disk model 
of Shakura & Sunvaev (1 19731) (^ = 4/3) and with the model of 
Agol & KroUk (.2000) (^ = 8/7), where the disk is powered by 
the spin of the centra l black hole. O ur result is less compati- 
ble with the model of iGaskell (|2008[) . which predicts a steeper 
energy profi le (4" - 1.75). Our resu lt compares well with that 
obtained by iPoindexter et al.l (l2008l) from multi-band photom- 
etry of the lensed quasar HE 1104-1805. They analyzed mi- 
crolensing lightcurves in eleven bands from the optical to the 
mid-infrared, and found that both the size and energy profile 
^ = 1 ■64^0 55 (i.e. their jS"') of the quasa r accretion disk are 
consistent with the standard disk model of I Shakura & SunvaevI 
(1I973.) . 

We observe that the determination of the power law-index is 
almost independent of the velocity prior used, whereas the de- 
termination of the scaled source size R, is directly proportional 
to the scaled effective transverse velocity V. This is easily ex- 
plained by the fact that the scaling between Rs and Rs is indepen- 
dent of the considered wavelength, and hence Rj/Rref = Rj/Rref- 
As a consequence, the determination of relative source sizes, and 
hence of the energy profile, is not expected to depend on the as- 
sumed microlens mass (M). 

Quasar microlensing is hence able to resolve structures of 
accretion disks on scales reaching 0.1 micro-arcsecond, which 
is more than 10, 000 times better than the resolution reached by 
today's best telescopes. Finally, we should mention several re- 
cent studies that compare microlensing in the X-ray and opti- 
cal domain, and that further demonstrate the efficiency of mi- 
crolensing in probing the inner parts of quasars. For instance, 
Morgan et al. (2008) analyze the lightcurves of the lensed quasar 
PG 1115+080, and find that the effective radius of the X-ray 
emission is 1.3^,j 5 dex smaller than that of the optical emis- 
sion, with the X-ray emission generated near the inner edge 
of the accretion disk while the optical emission comes from 
scales sli g htly lar ger than those expected for a standard thin 
disk. lPoolev et al.l (12008) observe extreme microlensing-induced 
(de-) magnification of the lensed images of PG 1115+080, and 
conclude that about 90% of the matter in the lensing galaxy 
is in smoothly distributed (dark) material and only about 10% 
is in compac t (stell ar) objects. Another example is given by 
Ch artas et al.l (120081) . who combine X-ray and optical data of 
HE 1104-1805, and reveal that the X-ray emitting region is 
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compact with a half-light radius smaller than six gravitational 
radius, i.e. smaller than 2 x lO'^ cm, thus placing significant 
constraints on active galactic nuclei models. 
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Table 6. Continuum integrated AB magnitudes in the six photometric bands defined in Table[2] 



ID HJD band 1 (mag) band 2 (mag) band 3 (mag) band 4 (mag) band 5 (mag) band 6 (mag) 



Image A 



1 


3292 


17.97 


± 


0.05 


17.81 


± 


0.05 


17.67 


± 


0.05 


17.54 


± 


0.06 


17.43 


± 


0.06 


17.33 


± 


0.05 


2 


3324 


18.02 


± 


0.05 


17.87 


± 


0.04 


17.73 


± 


0.04 


17.61 


± 


0.05 


17.50 


± 


0.05 


17.40 


± 


0.04 


3 


3341 


18.02 


± 


0.06 


17.85 


± 


0.05 


17.71 


± 


0.05 


17.59 


± 


0.06 


17.47 


± 


0.06 


17.37 


± 


0.05 


4 


3355 


17.98 


± 


0.08 


17.83 


± 


0.08 


17.69 


± 


0.07 


17.57 


± 


0.08 


17.47 


± 


0.08 


17.37 


± 


0.08 


5 


3502 


18.07 


± 


0.04 


17.88 


± 


0.04 


17.72 


± 


0.04 


17.57 


± 


0.04 


17.44 


± 


0.04 


17.33 


± 


0.04 


6 


3523 


18.16 


± 


0.05 


17.97 


± 


0.05 


17.80 


± 


0.04 


17.64 


± 


0.05 


17.51 


± 


0.05 


17.39 


± 


0.04 


7 


3553 


18.23 


± 


0.04 


18.03 


± 


0.04 


17.85 


± 


0.03 


17.69 


± 


0.04 


17.55 


± 


0.04 


17.42 


± 


0.03 


8 


3566 


18.16 


± 


0.06 


17.96 


± 


0.06 


17.79 


± 


0.05 


17.63 


± 


0.06 


17.49 


± 


0.06 


17.37 


± 


0.05 


9 


3589 


18.24 


± 


0.04 


18.04 


± 


0.04 


17.87 


± 


0.03 


17.71 


± 


0.04 


17.57 


± 


0.04 


17.45 


± 


0.03 


10 


3598 


18.24 


± 


0.06 


18.03 


± 


0.06 


17.84 


± 


0.06 


17.68 


± 


0.06 


17.53 


± 


0.06 


17.39 


± 


0.06 


11 


3608 


18.21 


± 


0.03 


18.02 


± 


0.03 


17.86 


± 


0.03 


17.71 


± 


0.04 


17.58 


± 


0.04 


17.47 


± 


0.03 


12 


3626 


18.21 


± 


0.05 


18.01 


± 


0.05 


17.83 


± 


0.04 


17.68 


± 


0.05 


17.54 


± 


0.05 


17.42 


± 


0.04 


13 


3641 


18.11 


± 


0.05 


17.90 


± 


0.05 


17.73 


± 


0.05 


17.57 


± 


0.05 


17.43 


± 


0.05 


17.30 


± 


0.05 


14 


3645 


18.08 


± 


0.05 


17.91 


± 


0.05 


17.77 


± 


0.05 


17.64 


± 


0.05 


17.52 


± 


0.05 


17.42 


± 


0.05 


15 


3655 


18.10 


± 


0.04 


17.92 


± 


0.04 


17.77 


± 


0.03 


17.63 


± 


0.04 


17.51 


± 


0.04 


17.39 


± 


0.03 


16 


3665 


18.10 


± 


0.04 


17.92 


± 


0.04 


17.75 


± 


0.04 


17.61 


± 


0.04 


17.48 


± 


0.04 


17.36 


± 


0.04 


17 


3686 


17.93 


± 


0.06 


17.77 


± 


0.06 


17.63 


± 


0.06 


17.50 


± 


0.07 


17.39 


± 


0.07 


17.28 


± 


0.06 


18 


3699 


17.94 


± 


0.04 


17.78 


± 


0.04 


17.64 


± 


0.04 


17.51 


± 


0.04 


17.40 


± 


0.04 


17.30 


± 


0.04 


19 


371 1 


17.88 


± 


0.07 


17.71 


± 


0.07 


17.57 


± 


0.07 


17.44 


± 


0.07 


17.33 


± 


0.07 


17.22 


± 


0.07 


20 


3880 


17.58 


± 


0.04 


17.47 


± 


0.05 


17.36 


± 


0.04 


17.26 


± 


0.05 


17.17 


± 


0.05 


17.08 


± 


0.04 


21 


3903 


17.50 


± 


0.03 


17.40 


± 


0.04 


17.30 


± 


0.04 


17.20 


± 


0.04 


17.12 


± 


0.04 


17.04 


± 


0.04 


22 


3907 


17.50 


± 


0.03 


17.39 


± 


0.03 


17.28 


± 


0.03 


17.19 


± 


0.04 


17.10 


± 


0.04 


17.02 


± 


0.03 


23 


3914 


17.51 


± 


0.02 


17.41 


± 


0.03 


17.30 


± 


0.03 


17.21 


± 


0.03 


17.12 


± 


0.03 


17.04 


± 


0.03 


24 


3944 


17.53 


± 


0.04 


17.39 


± 


0.04 


17.27 


± 


0.04 


17.16 


± 


0.04 


17.07 


± 


0.04 


16.98 


± 


0.04 


25 


3951 


17.51 


± 


0.03 


17.39 


± 


0.03 


17.28 


± 


0.03 


17.18 


± 


0.03 


17.10 


± 


0.04 


17.02 


± 


0.03 


26 


4022 


17.54 


± 


0.03 


17.41 


± 


0.04 


17.29 


± 


0.03 


17.19 


± 


0.04 


17.09 


± 


0.04 


17.01 


± 


0.03 


27 


4037 


17.53 


± 


0.03 


17.40 


± 


0.03 


17.29 


± 


0.03 


17.19 


± 


0.04 


17.09 


± 


0.04 


17.01 


± 


0.03 


28 


4050 


17.44 


± 


0.04 


17.33 


± 


0.04 


17.23 


± 


0.04 


17.14 


± 


0.04 


17.06 


± 


0.04 


16.99 


± 


0.04 


29 


4067 


17.41 


± 


0.05 


17.29 


± 


0.05 


17.19 


± 


0.05 


17.09 


± 


0.05 


17.01 


± 


0.05 


16.93 


± 


0.05 


31 


4092 


17.40 


± 


0.04 


17.30 


± 


0.05 


17.20 


± 


0.04 


17.12 


± 


0.05 


17.05 


± 


0.05 


16.99 


± 


0.04 


32 


4292 


17.51 


± 


0.03 


17.36 


± 


0.03 


17.22 


± 


0.03 


17.10 


± 


0.03 


16.99 


± 


0.03 


16.89 


± 


0.03 


33 


4297 


17.49 


± 


0.03 


17.33 


± 


0.03 


17.19 


± 


0.03 


17.06 


± 


0.03 


16.95 


± 


0.03 


16.85 


± 


0.03 


34 


4307 


17.42 


± 


0.03 


17.31 


± 


0.03 


17.21 


± 


0.03 


17.12 


± 


0.04 


17.04 


± 


0.04 


16.97 


± 


0.03 


35 


4316 


17.46 


± 


0.03 


17.34 


± 


0.03 


17.23 


± 


0.03 


17.13 


± 


0.03 


17.04 


± 


0.03 


16.97 


± 


0.03 


36 


4340 


17.43 


± 


0.04 


17.29 


± 


0.04 


17.17 


± 


0.04 


17.06 


± 


0.04 


16.97 


± 


0.04 


16.88 


± 


0.04 


37 


4350 


17.44 


± 


0.03 


17.36 


± 


0.03 


17.29 


± 


0.03 


17.23 


± 


0.03 


17.18 


± 


0.03 


17.13 


± 


0.03 


38 


4364 


17.40 


± 


0.03 


17.28 


± 


0.03 


17.18 


± 


0.03 


17.09 


± 


0.03 


17.00 


± 


0.03 


16.93 


± 


0.03 


39 


4367 


17.35 


± 


0.06 


17.23 


± 


0.06 


17.11 


± 


0.06 


17.02 


± 


0.07 


16.93 


± 


0.07 


16.85 


± 


0.06 


40 


4379 


17.39 


± 


0.02 


17.26 


± 


0.02 


17.15 


± 


0.02 


17.05 


± 


0.03 


16.96 


± 


0.03 


16.87 


± 


0.02 


41 


4384 


17.41 


± 


0.03 


17.29 


± 


0.03 


17.20 


± 


0.03 


17.11 


± 


0.04 


17.03 


± 


0.04 


16.96 


± 


0.03 


42 


4420 


17.39 


± 


0.05 


17.25 


± 


0.05 


17.13 




0.05 


17.02 


± 


0.05 


16.93 


± 


0.05 


16.84 


± 


0.05 


43 


4436 


17.39 


± 


0.03 


17.23 


± 


0.04 


17.09 


± 


0.03 


16.96 


± 


0.04 


16.85 


± 


0.04 


16.75 


± 


0.03 


Image B 


1 


3293 


19.08 


± 


0.12 


18.94 


± 


0.12 


18.82 


± 


0.11 


18.71 


± 


0.13 


18.62 


± 


0.13 


18.53 


± 


0.11 


2 


3324 


19.29 


± 


0.10 


19.11 


± 


0.09 


18.94 


± 


0.08 


18.79 


± 


0.09 


18.66 


± 


0.09 


18.54 


± 


0.09 


3 


3341 


19.18 


± 


0.14 


18.99 


± 


0.14 


18.83 


± 


0.12 


18.69 


± 


0.14 


18.57 


± 


0.14 


18.45 


± 


0.13 


4 


3356 


19.17 


± 


0.14 


19.00 


± 


0.14 


18.85 


± 


0.13 


18.72 


± 


0.14 


18.60 


± 


0.14 


18.49 


± 


0.13 


5 


3503 


18.43 


± 


0.05 


18.35 


± 


0.05 


18.28 


± 


0.04 


18.22 


± 


0.05 


18.17 


± 


0.05 


18.12 


± 


0.05 


6 


3523 


18.46 


± 


0.06 


18.36 


± 


0.06 


18.28 


± 


0.06 


18.21 


± 


0.06 


18.15 


± 


0.06 


18.09 


± 


0.06 


9 


3589 


18.76 


± 


0.05 


18.61 


± 


0.06 


18.48 


± 


0.05 


18.37 


± 


0.06 


18.26 


± 


0.06 


18.17 


± 


0.05 


11 


3608 


18.85 


± 


0.05 


18.68 


± 


0.05 


18.53 


± 


0.05 


18.40 


± 


0.05 


18.28 


± 


0.05 


18.18 


± 


0.05 


12 


3626 


18.67 


± 


0.06 


18.55 


± 


0.06 


18.44 


± 


0.06 


18.35 


± 


0.07 


18.26 


± 


0.07 


18.18 


± 


0.06 


13 


3641 


18.61 


± 


0.06 


18.48 


± 


0.07 


18.38 


± 


0.06 


18.28 


± 


0.07 


18.20 


± 


0.07 


18.12 


± 


0.06 


14 


3645 


18.62 


± 


0.07 


18.49 


± 


0.08 


18.38 


± 


0.07 


18.28 


± 


0.08 


18.19 


± 


0.08 


18.11 


± 


0.07 


15 


3655 


18.73 


± 


0.06 


18.58 


± 


0.06 


18.45 


± 


0.06 


18.34 


± 


0.06 


18.24 


± 


0.06 


18.15 


± 


0.06 


16 


3665 


18.72 


± 


0.06 


18.58 


± 


0.06 


18.46 


± 


0.06 


18.35 


± 


0.07 


18.25 


± 


0.07 


18.16 


± 


0.06 


17 


3686 


18.63 


± 


0.10 


18.50 


± 


0.10 


18.38 


± 


0.09 


18.28 


± 


0.10 


18.19 


± 


0.10 


18.10 


± 


0.09 


18 


3699 


18.66 


± 


0.08 


18.53 


± 


0.08 


18.41 


± 


0.08 


18.31 


± 


0.09 


18.22 


± 


0.09 


18.13 


± 


0.08 


19 


3711 


18.56 


± 


0.11 


18.42 


± 


0.11 


18.30 


± 


0.10 


18.20 


± 


0.11 


18.10 


± 


0.11 


18.01 


± 


0.10 


20 


3880 


18.83 


± 


0.12 


18.74 


± 


0.13 


18.59 


± 


0.12 


18.47 


± 


0.14 


18.35 


± 


0.14 


18.25 


± 


0.12 


21 


3903 


18.89 


± 


0.07 


18.77 


± 


0.07 


18.60 


± 


0.07 


18.45 


± 


0.08 


18.31 


± 


0.08 


18.19 


± 


0.07 


22 


3907 


18.83 


± 


0.06 


18.71 


± 


0.07 


18.53 


± 


0.07 


18.38 


± 


0.07 


18.25 


± 


0.07 


18.12 


± 


0.07 


23 


3914 


18.86 


± 


0.05 


18.74 


± 


0.06 


18.57 


± 


0.05 


18.42 


± 


0.06 


18.29 


± 


0.06 


18.17 


± 


0.05 


24 


3944 


18.76 


± 


0.07 


18.60 


± 


0.07 


18.47 


± 


0.06 


18.35 


± 


0.07 


18.24 


± 


0.07 


18.15 


± 


0.06 


25 


3951 


18.72 


± 


0.06 


18.56 


± 


0.06 


18.43 


± 


0.06 


18.31 


± 


0.06 


18.20 


± 


0.06 


18.10 


± 


0.06 


26 


4022 


18.76 


± 


0.06 


18.61 


± 


0.06 


18.48 


± 


0.05 


18.36 


± 


0.06 


18.25 


± 


0.06 


18.15 


± 


0.05 


27 


4037 


18.69 


± 


0.05 


18.54 


± 


0.05 


18.41 


± 


0.05 


18.29 


± 


0.05 


18.19 


± 


0.05 


18.09 


± 


0.05 


28 


4050 


18.63 


± 


0.09 


18.49 


± 


0.09 


18.36 


+ 


0.08 


18.25 


± 


0.09 


18.15 


± 


0.09 


18.06 


± 


0.08 


29 


4067 


18.59 


± 


0.12 


18.46 


± 


0.12 


18.34 


+ 


0.12 


18.23 


± 


0.12 


18.13 


± 


0.12 


18.05 


± 


0.12 


30 


4089 


18.57 


± 


0.14 


18.45 


± 


0.14 


18.34 


± 


0.13 


18.25 


± 


0.14 


18.16 


± 


0.14 


18.08 


± 


0.13 


31 


4093 


18.60 


± 


0.10 


18.46 


± 


0.10 


18.34 


± 


0.09 


18.23 


± 


0.10 


18.13 


± 


0.10 


18.05 


± 


0.09 


32 


4292 


18.66 


± 


0.06 


18.49 


± 


0.07 


18.34 


± 


0.06 


18.21 


± 


0.07 


18.09 


± 


0.06 


17.98 


± 


0.06 


33 


4297 


18.66 


± 


0.07 


18.49 


± 


0.07 


18.35 


± 


0.06 


18.22 


± 


0.07 


18.11 


± 


0.07 


18.00 


± 


0.06 


34 


4307 


18.54 


± 


0.09 


18.40 


± 


0.09 


18.28 


± 


0.08 


18.17 


± 


0.09 


18.08 


± 


0.09 


17.99 


± 


0.08 


35 


4316 


18.66 


± 


0.08 


18.50 


± 


0.09 


18.36 


± 


0.08 


18.24 


± 


0.09 


18.12 


± 


0.09 


18.02 


± 


0.08 


36 


4340 


18.69 


± 


0.13 


18.51 


± 


0.13 


18.35 


± 


0.12 


18.20 


± 


0.13 


18.07 


± 


0.13 


17.96 


± 


0.12 


37 


4350 


18.59 


± 


0.08 


18.46 


± 


0.08 


18.34 


± 


0.08 


18.23 


± 


0.08 


18.14 


± 


0.08 


18.05 


± 


0.08 


38 


4364 


18.60 


± 


0.09 


18.44 


± 


0.09 


18.31 


± 


0.09 


18.20 


± 


0.10 


18.09 


± 


0.09 


18.00 


± 


0.09 


39 


4367 


18.59 


± 


0.16 


18.46 


± 


0.17 


18.34 


± 


0.15 


18.24 


± 


0.17 


18.15 


± 


0.17 


18.07 


± 


0.16 


40 


4379 


18.62 


± 


0.05 


18.45 


± 


0.06 


18.30 


± 


0.05 


18.17 


± 


0.06 


18.05 


± 


0.06 


17.94 


± 


0.05 


41 


4384 


18.68 


± 


0.06 


18.51 


± 


0.07 


18.36 


± 


0.06 


18.24 


± 


0.07 


18.12 


± 


0.07 


18.01 


± 


0.06 


42 


4421 


18.72 


± 


0.09 


18.54 


± 


0.09 


18.39 


± 


0.08 


18.25 


± 


0.09 


18.13 


± 


0.09 


18.02 


± 


0.08 


43 


4436 


18.60 


± 


0.07 


18.43 


± 


0.07 


18.28 


± 


0.06 


18.14 


± 


0.07 


18.02 


± 


0.07 


17.91 


± 


0.06 
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Table 6. Continued. 



ID 


HJD 


band 1 (mag) 


band 2 (mag) 


band 3 (mag) 


band 4 (mag) 


band 5 (mag) 


band 6 (mag) 


Image C 




jzyj 


1 o no ^ n 1 1 


1 fi 07 J. n 17 

icS.yz ± u. iz 


1 Q 713 ^ n 1 1 
Lo./o ± u.l I 


1 8 ^ n 1 1 

Io.Oj ± U.l J 


1 8 ^/i J. n 1 1 
lo.j4 ± u.l J 


1 9 AA J. n 1 1 
15.44 ± U.ll 


2 


^JjZ'f 


1 Q "20 J. n no 
iy.:>z ± u.uy 


1 Q no -u n nQ 


1 Q on 4. n ns 

lo.VU ± U.Uo 


1 8 71 -1- n no 
lo./o ± u.u? 


1 8 S7 -1- n no 

lo. J / ± u.uy 


1 9 Al -1- n na 

15.4J ± U.U5 


-1 

D 




1 Q 70 -L n 1 /I 

ly.zz ± u. 14 


1 Q A1 -L A 1/1 

ly.ui ± u. 14 


1 Q Q7 -1- A 1 1 
1 fi.fiZ ± U. I J 


1 Q -1- n 1 /I 

1 O.OJ ± U. 14 


1 Q '^n -L n 1 /I 

Io.jU ± U. 14 


1 9 17 -1- A 1 7 
lo.J / ± V/. IZ 


4 


3356 


1 Q AA n 1 A 
1 y.44 ± U. 1 D 


1 77 _i_ A 1 A 
ly.zz ± U. lO 


1 n7 -1- A 1 A 
IV.UZ ± u. lo 


1 Q sj^ -1- n 1 ft 

1 o.o J ± u.lo 


1 ft AO n 1 ft 
1 o.oy ± u.lo 


1ft «;s A IS 

lo.JJ ± \). 1 J 


5 


3503 




1 Q AQ -t- A AX 


1 X 01 + A AS 


1 Si 70 -1- n no 

i o. / y i u.uy 


1 ft AA _(_ n nft 

1 O.UO i V/.V/O 


1ft S4 + A AS^ 

lO.JH- i U.UO 




3523 


1 7Q _L n 1 A 

1 v.zv ± u. ivy 


1 AQ -L A 1 n 
1 7.\}y ± \). iU 


1 Q 1 -1- A 1 n 
1 o.y 1 ± u. lu 


1 Si 7A -1- n 11 
1 o. / ± U. 1 1 


1 ft A7 n 1 A 
1 o.OZ ± V). L\) 


1ft ^n A AO 

1 o. jvj ± U.uy 


Q 


JJOV 


1 Q A^L -t. n no 
iy.40 ± u.uy 


10 71 J. n no 
ly.zj ± u.uy 


1 n7 -1- n nft 
ly.uz ± U.Uo 


1 8 81 -1- n no 
lo.oj ± u.uy 


1 s A7 -I. n no 
lo.D/ ± u.uy 


19 S7 -•- n na 

lo.JZ ± U.Uo 


1 1 


OUUo 


1 Q ^0 J. n no 


1 Q 19 J. n no 
ly. jz ± u.uy 


10 1 n 4. n ns 
ly. lu ± U.Uo 


1 8 8Q 4. n no 
15. 5y ± u.uy 


18 71 -1- n ns 

lo. /I ± U.Uo 


1 9 SA -1- n na 

15.J4 ± U.U5 


12 


3626 




19 17 ± 11 


1 R QQ ^. n in 
Lo.yy I u. lu 


IR Ri J. n 11 

16. 6 J i U. 1 1 


1 S fiR + n 1 1 

10.06 I U. 1 1 


1 R ss + n in 

10. X U. lU 


13 


3641 


1Q 91 + n in 


1 Q n9 + n 1 n 

l3.UZ I U. lU 


1 R R^ J. n in 

lo.oJ i U. lU 


1 R 7n -1- n 11 

16. /U i U. 1 1 


1 s sfi + n 1 

lo.JD X U. lU 


1 R /M + n nQ 
10.44 X U.uy 


14 




1 Q 0^ J. n 1 

ly.zj ± u. iz 


1 Q f\A -1- n 1 9 
13.U4 ± U. IZ 


IS 4. n 19 

lo.oO ± U. IZ 


1 8 7n J. n 1 1 
15. /u ± u. lo 


1 8 SA -1- n 1 1 

lo.DO ± yj.LD 


1 a Al -1- n 1 1 

16.4J X U. 1 1 


15 


3655 




1 90 + A in 


1 Q nn + A AQ 


1 Q Q9 + n 10 

1 o.oz i u. IV/ 


1 ft A< 4_ n in 

1 O.U J i V/. IV/ 


19 sn n nQ 

lo.jv/ u.uy 


iO 




1 y.4 1 ± U. 


10 1 -(- A in 
1 y. 1 V ± u. 1 u 


1 « 00 -1- A AQ 

1 o.yy ± u.uv 


1 Si Si 1 -1- n 1 A 
1 o.o 1 ± U. l\) 


1 ft AA n 1 A 
1 o.Oo ± V). L\) 


1ft ^1 A AO 

1 O.J 1 ± U.uy 


17 


3686 


1 y. i V ± u. 1 J 


1 QSl _t_ A 1 s 

lo.yo ± u. 1 J 


1 « Si A -1- A 1 zl 
1 o.oU ± U. 1 4 


1 Si Al -1- n 1 s 

1 o.OJ ± U. 1 J 


1 ft Asi n 1 ^ 

1 o.4o ± VJ. 1 J 


1 ft IS A 1 1 

lo.JJ ±\).LD 


16 


3699 


1 74 n 17 

1 y.z4 ± u. 1 z 


1 Al A 17 

1 y.uj ± u. 1 z 


1 SI SIA -1- A 1 1 
1 o.o4 ± U. 1 1 


1 Si Aft -1- n 17 
1 o.Oo ± U. 1 Z 


1ft -L n 17 

1 O.J J ± VJ. 1 Z 


1 ft An -L A 1 1 

1 0.4U ± VA 1 I 


1 Q 


^7 1 1 


1 fi OA 4- n 1 zi 

Lo.y^ ± U. 14 


1 fi lA -t. A 1 A 
lo. /4 ± U. 14 


1 Q <7 a- A 1 A 
1 0. J / ± U. 14 


1 8 A7 J- n 1 <^ 
1 0.4Z ± U. Lj 


1 S 78 -1- n 1 <^ 
lO.Zo ± U. 1 J 


1ft 1 A -i- A 11 
lo.lO ± U.l J 




JOOU 


lo.oO ± U.lo 


1 Q ifi J. n 9n 

lo. /O ± u.zu 


1 8 /;i -1- n 9n 

15.01 ± U.ZU 


1 8 A7 J. n 99 
15.4/ ± U.ZZ 


1 a IS -1- n 91 

lO.JJ ± U.Zl 


1 a 9A -1- n 1 
10. Z4 X u. ly 


21 




ly. ly ± u. iz 


1 Q n^ -1- n 1 A 
ly.uj ± u. 14 


18 8A 4. n 11 
15.54 ± U.l J 


1 8 ft7 -1- n 1 A 
15.0/ ± U.14 


1 a SI -1- n 1 A 

lo.Jl ± u.14 


1 a 1A -1- n 1 1 

10.30 X U.lJ 


22 


3907 


1Q 19 + n 1 

13. IZ I u. u 


1 Q OR + n 14 

10.30 X U. 14 


1 R 77 J. n 14 

16. / / ± U. 14 


1 R SO -1- n 1 s 

16. jy ± U. 1 J 


1 R ill + n 1 s 

10.4J X U. 1 J 


1 R 9Q + n 11 
10. zy X u. u 






ly.ij ± u.iu 


1 Q ni J. n 1 7 
ly.ui ± u.iz 


IB 81 j.n 11 

lo.ol ± U. 1 1 


1 8 #;i J. n n 

lo.OO ± U.IZ 


1 a y17 J. n 17 
lo.4/ ± U. IZ 


1 a 17 J. n 1 1 
10. oz ± u.l 1 




"20/1/1 


1 Q ni -1- n CiQ 

IV. ul ± U.Uo 


1 Q Q'i -1- A AO 


1 Q A7 -1- A AQ 
15. D / ± U.Uo 


I Si =^7 -1- n no 
1 fs. jz + U.uy 


1 9 10 J. n no 
lo.jy ± u.uy 


19 7 Si -1- A ASi 
lo.ZO ± \J.\JIS 


Zj 


JVj 1 


1 o.yv ± u.u / 


1 ft ftn A ASi 
lo.oU ± U.Uo 


1 Q AA -1- A A7 
1 0.04 ± U.U / 


1 Si ^A -1- n nft 

1 o.ZiU ± U.Uo 


1 ft 17 _L n n7 

lo.J / ± \.).\.) 1 


1 ft -(- A A7 
lo.ZJ ± u.u / 


ZD 


4022 


10 1 zi n 11 
1 y. 1 4 ± u. 1 1 


1 ft Ozl A 1 7 

io.y4 ± u. iz 


1 9 7A -1- A 1 1 
1 o. / ± U. 1 1 


1 Si Ai -1- n 1 7 
lo.ol ± u. iz 


1 S AA n 1 7 
lo.4D ± U.IZ 


1 ft 1A -1- A in 
lo.J4 ± U. iK) 


27 


4037 


1 9 07 ± 11 


1 S S7 -t- A 1 1 


1 X 7A + A 1 n 

1 O. /V/ X V/. 1 


1 Si ss + n 1 1 

1 O.J J i vj. i 1 


18 41 ± 1 1 


1 ft 7Q _(_ A in 

I O.ZV V/. 1 V/ 


ZO 




1 Q SI ..i. n 1 1 
lo.oJ ± U. 1 1 


1 S -i. A 1 1 
lo.DD ± U.l 1 


1 Q <f\ J.. A 1 1 
lO.DU ± U.l 1 


1 Si 1A J. n 1 7 
lo.JD ± U.IZ 


19 7/1 J., n 11 

15. Z4 ± U. 1 1 


19 1 7 a, A 11 
lo. IZ ± U.l 1 


zy 


'+UO/ 


lo. /o ± u.iy 


1 Q Al -1- n 7n 

lo.Ol ± u.zu 


1 8 A7 4. n 1 
15.4/ ± U. l7 


1 8 1A -1- n 9n 

15. J4 ± U.ZU 


1 a 99 -1- n 9n 

15.ZZ ± u.zu 


19 11 -unio 
10. 1 1 X u. ly 




'+U07 


1 Q fin -1- n 91 

lo.oU ± U.Zl 


1 S A9 -1- n 91 
lo.OZ ± U.ZJ 


1 8 A7 -1- n 99 
15.4 / ± U.ZZ 


1 8 11 -1- n 91 

15. J J ± U.ZO 


18 91 -1- n 91 
15. Zl ± U.ZO 


19 1 n -1- n 9n 
10. lu X u.zu 


1 




1 Q Q9 -1- n 9n 
lo.yz ± u.zu 


1 9 79 -1- n 9n 

lo. /Z ± u.zu 


18 -1. n 10 

lo.JJ ± U.iy 


1 8 An -1- n 9 1 

1 0.4U ± U.Zl 


1 9 77 -1- n 9n 

lo.Z/ ± u.zu 


19 1S-i-ni9 
10.1 J X U.16 


32 


4292 


18.75 ± 0.06 


18.57 ± 0.06 


18.42 ±0.06 


18.28 ±0.06 


18.16 ± 0.06 


18.05 ± 0.06 


33 


4297 


18.73 ±0.06 


18.57 ±0.07 


18.42 ±0.06 


18.29 ±0.07 


18.18 ±0.07 


18.07 ±0.06 


34 


4307 


18.57 ±0.08 


18.43 ±0.09 


18.31 ±0.08 


18.20 ±0.09 


18.11 ±0.09 


18.02 ±0.08 


35 


4316 


18.71 ±0.08 


18.55 ±0.09 


18.42 ±0.08 


18.29 + 0.09 


18.18 ±0.09 


18.08 ±0.08 


36 


4340 


18.79 ±0.12 


18.59 + 0.13 


18.42 ± 0.12 


18.26 + 0.13 


18.12 ±0.13 


17.99 + 0.12 


37 


4350 


18.79 ± 0.09 


18.63 ±0.09 


18.50 ±0.08 


18.38 + 0.09 


18.27 ± 0.09 


18.18 ±0.08 


38 


4364 


18.80 ± 0.09 


18.62 ±0.10 


18.46 ± 0.09 


18.32 ±0.10 


18.19 ±0.10 


18.08 ±0.09 


39 


4367 


18.50 ±0.14 


18.37 ±0.15 


18.26 ±0.14 


18.17 ±0.16 


18.08 ± 0.16 


18.00 ±0.14 


40 


4379 


18.98 ±0.07 


18.77 ±0.08 


18.58 + 0.07 


18.41 +0.08 


18.26 ±0.07 


18. 12 ±0.07 


41 


4384 


19.03 ±0.09 


18.82 ±0.10 


18.64 ± 0.09 


18.48 + 0.10 


18.34 ±0.10 


18.21 ±0.09 


42 


4421 


19.18 ±0.13 


18.94 ±0.13 


18.74 ±0.12 


18.55 ±0.14 


18.39 ±0.13 


18.24 ±0.12 


43 


4436 


19.06 ± 0.09 


18.83 ±0.10 


18.63 ±0.10 


18.45 ±0.10 


18.29 ±0.10 


18.14 ±0.09 



Image D 



1 


3292 


19.11 ±0.14 


18.91 ±0.14 


18.73+0.14 


18.57 + 0.15 


18.42 ±0.15 


18.29 ±0.13 


2 


3324 


19.23 ±0.12 


19.03 ±0.12 


18.86 ±0.12 


18.71 ±0.13 


18.58 ±0.13 


18.45 ±0.12 


3 


3341 


19.18 ±0.14 


18.98 ±0.14 


18.79 ±0.14 


18.63 ±0.15 


18.49 ±0.15 


18.36 ±0.13 


4 


3355 


19.04 ±0.18 


18.84 ±0.19 


18.67 ±0.18 


18.51 ±0.21 


18.38 ± 0.20 


18.25 ±0.18 


5 


3502 


19.25 ±0.11 


19.02 ±0.11 


18.81 ±0.10 


18.63 ±0.11 


18.47 ±0.11 


18.32 ±0.10 


6 


3523 


19.29 ±0.11 


19.08 ±0.11 


18.89 ±0.11 


18.73 ±0.12 


18.58 ± 0.12 


18.45 ±0.11 


7 


3553 


19.42 ±0.10 


19.20 ±0.10 


19.01 ±0.09 


18.84 ±0.11 


18.68 ±0.10 


18.55 ±0.09 


8 


3566 


19.26 ±0.14 


19.01 ±0.14 


18.80 ±0.13 


18.61 ±0.15 


18.44 ±0.14 


18.28 ±0.13 


9 


3589 


19.44 ±0.09 


19.23 ±0.09 


19.04 ± 0.09 


18.87 ±0.10 


18.73 ±0.10 


18.59 ±0.09 


10 


3598 


19.35 ±0.15 


19.10 ±0.15 


18.87 ±0.14 


18.68 ±0.15 


18.50 ±0.15 


18.34 ±0.13 


11 


3608 


19.43 ± 0.09 


19.22 ± 0.09 


19.05 ± 0.08 


18.89 ±0.10 


18.75 ±0.10 


18.62 ± 0.09 


12 


3626 


19.43 ±0.12 


19.20 ±0.12 


18.99 ±0.11 


18.81 ±0.12 


18.65 ±0.12 


18.50 ±0.11 


13 


3641 


19.27 ±0.12 


19.02 ±0.12 


18.79 ±0.11 


18.60 ±0.13 


18.42 ±0.13 


18.26 ±0.11 


14 


3645 


19.26 ±0.12 


19.03 ±0.12 


18.83+0.11 


18.65 ±0.13 


18.50 ±0.13 


18.35 ±0.11 


15 


3655 


19.33 ±0.10 


19.13 ±0.10 


18.96 + 0.09 


18.81 ±0.11 


18.67 ±0.11 


18.54 ±0.10 


16 


3665 


19.33 ±0.11 


19. 10 ±0.11 


18.91 +0.10 


18.73 ±0.12 


18.58 ±0.11 


18.43 ±0.10 


17 


3686 


19.03 ±0.15 


18.82 ±0.15 


18.64 ±0.15 


18.48 ±0.16 


18.33 ±0.16 


18.20 ±0.14 


18 


3699 


19.14±0.11 


18.94 ±0.11 


18.76±0.11 


18.60 ±0.12 


18.46 ±0.12 


18.33 ±0.10 


19 


3711 


18.95 ±0.17 


18.75 ±0.18 


18.58 ±0.17 


18.42 ±0.19 


18.28 ±0.18 


18.16 ±0.16 


20 


3880 


19.00 ±0.11 


18.83 ±0.13 


18.66 ±0.12 


18.51 ±0.14 


18.37 ±0.13 


18.25 ±0.12 


21 


3903 


19.08 ±0.10 


18.93 ±0.12 


18.77 ±0.11 


18.63 ±0.13 


18.50 ±0.13 


18.39 ±0.11 


22 


3907 


19.07 ± 0.08 


18.92 ± 0.09 


18.75 ±0.09 


18.61 ±0.10 


18.48 ± 0.10 


18.36 ± 0.09 


23 


3914 


19.13 ±0.07 


18.98 ±0.08 


18.83 ±0.07 


18.70 ±0.08 


18.58 ±0.08 


18.47 ±0.08 


24 


3944 


19.08 ±0.12 


18.87 ±0.12 


18.69 ± 0.12 


18.52 ±0.13 


18.38 ±0.13 


18.24 ±0.11 


25 


3951 


19.04 ±0.11 


18.85 ±0.10 


18.68 ± 0.10 


18.53 + 0.11 


18.40 ±0.11 


18.28 ±0.10 


26 


4022 


19.06 ±0.07 


18.88 ±0.08 


18.72 ± 0.07 


18.58 + 0.08 


18.46 ±0.08 


18.34 ±0.07 


27 


4037 


19.03 + 0.07 


18.85 + 0.07 


18.70 ± 0.07 


18.56 + 0.08 


18.44 + 0.08 


18.32 + 0.07 


28 


4050 


18.93 + 0.12 


18.75 + 0.12 


18.59 ±0.12 


18.45+0.13 


18.32 + 0.13 


18.20 + 0.12 


29 


4067 


18.84 ±0.11 


18.65 ±0.11 


18.49 ±0.10 


18.35 + 0.12 


18.22 + 0.12 


18.11 ±0.10 


31 


4092 


18.87 ± 0.09 


18.70 ± 0.09 


18.55 ± 0.09 


18.42 ±0.10 


18.30 ±0.10 


18.19 ±0.09 


32 


4292 


19.00 ± 0.09 


18.78 ± 0.09 


18.58 ±0.09 


18.41 ±0.10 


18.25 ± 0.09 


18.11 ±0.08 


33 


4297 


18.95 ±0.08 


18.73 ±0.08 


18.53 ±0.08 


18.36 ±0.09 


18.21 ±0.08 


18.07 ±0.08 


34 


4307 


18.81 ±0.11 


18.63 ±0.12 


18.48 ± 0.12 


18.34 ±0.13 


18.22 ±0.13 


18.11 ±0.11 


35 


4316 


18.85 ±0.11 


18.69 ±0.12 


18.54 ±0.11 


18.41 ±0.12 


18.30 ±0.12 


18.20 ±0.11 


36 


4340 


18.73 ±0.13 


18.55 ±0.13 


18.39 ±0.13 


18.24 ±0.14 


18.12 ±0.14 


18.00 ±0.12 


37 


4350 


18.86 ±0.11 


18.75 ±0.12 


18.65 ±0.12 


18.57 ±0.14 


18.50 ±0.14 


18.43 ±0.12 


38 


4364 


18.76 ±0.10 


18.59 ±0.11 


18.44 ±0.10 


18.31 ±0.11 


18.19 ±0.11 


18.09 ±0.10 


39 


4367 


18.59 ±0.18 


18.41 ±0.19 


18.26 ±0.18 


18.12 ±0.20 


17.99 ± 0.20 


17.88 ±0.18 


40 


4379 


18.70 ± 0.07 


18.52 ± 0.08 


18.36 ± 0.08 


18.22 ± 0.08 


18.09 ± 0.08 


17.97 ±0.07 


41 


4384 


18.74 ± 0.09 


18.57 ±0.10 


18.42 ±0.10 


18.29 ±0.11 


18.18 ±0.11 


18.07 ± 0.09 


42 


4420 


18.61 ±0.12 


18.43 ±0.13 


18.27 ±0.13 


18.13 ±0.14 


18.01 ±0.14 


17.89 ±0.13 


43 


4436 


18.59 ± 0.09 


18.39 ±0.10 


18.21 ±0.09 


18.05 ±0.10 


17.91 ±0.10 


17.78 ± 0.09 



